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A  FOLLY  MACROSCOPIC  DESCRIPTION  OF  BOUNDED  SEMICONDUCTORS  WITH 
«  AN  APPLICATION  TO  THE  Si-Si02  INTERFACE 

M.G.  Ancona  and  H.F.  Tiersten 
Department  of  Mechanical  Engineering, 

Aeronautical  Engineering  &  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 

ABSTRACT 

A  fully  macroscopic  description  of  semiconductors  is  presented 
which  includes  the  boundary  conditions  at  the  surface  of  the  semicon¬ 
ductor  that  are  required  for  consistency  with  the  usual  diffusion- 
drift  current  equations.  As  in  all  field  theories,  e.g.,  electro¬ 
magnetism,  both  the  boundary  conditions  and  the  differential  equations 
are  obtained  from  the  same  governing  integral  forms.  The  new  boundary 
conditions  relate  the  jump  discontinuities  in  the  chemical  potentials 
across  the  interface  to  the  forces  exerted  by  the  lattice  on  the  charge 
carriers  which  prevent  the  carriers  from  leaving  the  solid.  The  expres¬ 
sions  for  the  forces  in  the  static  case  are  found  and  the  values  of  the 
,  material  surface  coefficients  appearing  therein  are  obtained  from  quasi¬ 

static  MOS  c-V  measurements  for  some  particular  Si-Si02  interfaces. 


1.  Introduction 


The  basic  understanding  of  the  behavior  of  semiconductors  arises  from 

the  semiclassical  description  of  the  motion  of  electrons  in  a  periodic 
1-3 

potential  .  This  model  of  the  semiconductor  has  proven  invaluable  in 

interpreting  a  host  of  phenomena  in  the  infinite  lattice.  Nevertheless, 

4 

when  the  electron  number  density  varies  appreciably  ,  as  it  often  does  near 

boundary  surfaces,  or  in  cases  of  inhomogeneous  doping,  a  macroscopic 

description,  which  is  not  related  in  a  precise  quantitative  way  to  the 

5-7 

above-mentioned  microscopic  model,  is  employed  .  This  standard  macro- 

Q 

scopic  description  ,  consisting  of  diffusion-drift  current  equations,  the 

charge  balance  equations  and  the  electrostatic  constitutive  equations,  is  a 

continuum  field  theory,  i.e.,  a  system  of  partial  differential  equations 

containing  dependent  macroscopic  field  variables  as  a  function  of  space  and 

time.  However,  a  complete  field  theory  must  also  have  a  set  of  consistent 

boundary  conditions.  These  are  obtained  from  integral  forms  of  the  governing 

equations  just  as  in  electromagnetism,  where  electromagnetic  boundary  con- 

9-10 

ditions  are  obtained  from  the  integral  forms  of  Maxwell's  equations  .  In 

contrast  to  electromagnetism,  however,  the  integral  forms  and  the  associated 

boundary  conditions  have  never  been  obtained  for  the  standard  macroscopic 

5-7 

theory  of  semiconductors 

A  well-defined  macroscopic  model  of  a  semiconductor  has  been  presented 
in  the  literature11  along  with  the  integral  forms  obtained  from  that  model. 
The  integral  forms  have  been  shown  to  yield  the  conventional  diffusion- 
drift  current  differential  equations  plus  the  associated  boundary  conditions 
across  the  surface  of  the  semiconductor,  which  are  missing  in  the  existing 
macroscopic  semiconductor  theory.  However,  since  the  equations  appearing 
in  Ref. 11  are  unduly  complicated  by  the  inclusion  of  mechanical  deformations 


2 


in  the  description  and,  for  other  reasons,  are  not  in  a  convenient  form,  we 

rederive  the  equations  and  obtain  the  more  useful  form  here,  Hie  macroscopic 

12 

model  from  which  the  field  theory  follows  is  chosen  to  permit  electronic  and 
hole  conduction  in  interaction  with  the  electrostatic  field.  Although  the 
model  can  readily  be  extended  to  incorporate  electric  polarization1*,  we 
do  not  bother  with  this  refinement  and  include  polarization  in  the  usual 
linear  way.  Also,  in  the  interest  of  clarity  all  resonance  phenomena  are 
expressly  excluded  from  the  treatment  by  ignoring  the  mass  of  all  moving 
components  in  the  model.  Accordingly,  the  model  consists  of  three  suitably 
defined  interpenetrating  continua.  Which  are  called  the  lattice  continuum, 
the  conduction  electronic  continuum  and  the  hole  continuum.  Hie  lattice  con¬ 
tinuum  does  not  move  and  contains  a  charge  density  representative  of  the 
fixed  ionized  impurities  responsible  for  extrinsic  semiconduction.  The  con¬ 
duction  electronic  and  hole  continua  are  inertialess  charged  fluids  which 
move  through  the  lattice  continuum  while  experiencing  a  force  of  resistance. 

In  addition,  each  conducting  fluid  interacts  with  neighboring  elements  of 
the  same  fluid  by  means  of  defined  fluid  pressure  forces.  Furthermore,  charge 
exchange  is  permitted  between  all  three  continua  in  order  to  allow  for  re¬ 
combination-generation  phenomena. 

As  in  all  such  descriptions,  the  application  of  the  laws  of  balance  of 

charge  and  momentum  yields  the  equations  of  motion  of  the  matter,  which, 

with  the  equations  of  electrostatics,  constitutes  an  underdetermined  system. 

Hie  system  is  completed  by  the  addition  of  various  materially  descriptive 

constitutive  relations11  whose  allowed  functional  forms  are  determined  from 

14 

thermodynamic  argtnents  .  Hie  differential  forms  of  the  equations  of  motion 


for  the  conduction  electronic  and  hole  fluids  along  with  linear  constitutive 


assumptions  yields  precisely  the  usual  diffusion-drift  current  equations 
containing  the  material  coefficients  of  mobility  and  diffusivity. 


The  boundary  conditions  at  the  surface  of  the  semiconductor,  which  are 
required  for  the  solution  of  semiconduction  problems  involving  bounded  media, 
are  obtained  from  integral  forms  of  the  semiconduction  equations.  As  with  all  the 
other  differential  equations  and  their  associated  boundary  conditions,  this 
procedure  ensures  consistency  between  the  semiconduction  differential  equa¬ 
tions  and  the  semiconduction  boundary  conditions.  These  boundary  conditions, 

1-7 

which  do  not  appear  in  any  of  the  existing  literature  on  semiconductors  ’ 
relate  the  jump  discontinuities  in  the  electronic  and  hole  chemical  potentials 
across  the  surface  to  the  forces  per  unit  area  per  unit  charge  densities 
exerted  by  the  lattice  continuum  on  the  respective  fluids  that  keep  the 
electrons  from  leaving  the  solid15.  The  thermodynamic  arguments  in 
Appendix  A  indicate  that  the  forces  per  unit  area  per  unit  charge  densities 
are  functions  of  the  electric  field  on  each  side  of  the  surface  of  the  semi¬ 
conductor  and  the  respective  charge  densities  of  each  fluid  at  the  surface. 

The  expressions  for  the  forces  per  unit  area  per  unit  charge  densities  con¬ 
tain  material  surface  coefficients  which  are  to  be  determined  from  measure¬ 
ments  just  as  the  mobility  and  diffusivity  coefficients  in  the  differential 
equations  for  the  semiconductor  are  often  measured.  Presumably,  the  material 
surface  coefficients  could  be  calculated  from  a  more  fundamental- 

quantum  mechanical  model  by  means  of  electronic  surface  structure  calcu- 
16  17 

lations  ’  ,  but  the  available  results  seem  to  indicate  that  for  quanti- 

18  19 

tative  detail  such  calculations  would  be  prohibitively  complicated  *  . 


A 


4. 


The  missing  boundary  conditions  have  not  deterred  workers  from  treating 
problems  concerning  semiconductors  with  boundaries  analytically  .  To  be 
sure,  many  authors  have  circumvented  the  difficulty  by  imposing  various 
boundary  conditions  in  order  to  provide  the  additional  condition  required 
for  the  solution  of  the  boundary  value  problem.  For  example,  some  have 
simply  disregarded  the  existence  of  surface  charge  on  the  semiconductor 


surface 


21 


Others,  especially  those  interested  in  computer  aided 


22  23-25 

design  have  assumed  apriori  the  value  of  the  surface  charge  .  Another 

26  27 

attack  has  been  that  of  Kroemer  *  ,  who  used  the  concept  of  a  "control  char¬ 

acteristic"  in  treating  certain  aspects  of  Gunn  effect  phenomena.  Many 

28 

additional  examples  may  be  cited  .  Although  some  of  these  procedures  are 

29 

based  on  reasonable  assumptions  over  certain  ranges  ,  in  a  number  of  cases 

.  30 

results  have  been  obtained  which  are  clearly  at  variance  with  experiment 

In  addition,  none  of  the  previous  approaches  have  formulated  the  boundary 
value  problem  without  assuming  boundary  conditions  not  deduced  from  funda¬ 
mental  principles.  In  the  procedure  presented  here  when  the  newly  defined 
material  surface  coefficients  have  been  found  from  measurements  of  a  partic¬ 
ular  material  surface,  the  new  boundary  conditions  enable  the  solution  of 
semiconduction  boundary  value  problems. 

In  the  final  section  of  this  paper  an  illustration  of  this  purely 
macroscopic  description  of  surfaces  of  semiconductors  is  presented  by 
treating  the  often  analyzed  semiconductor- insulator  interface.  Much  work' 
has  been  done  with  the  purpose  of  understanding  and  describing  the 

semiconductor-insulator  junction,  especially  in  the  case  of  the  Si-Sio2 
34 

interface  .  Microscopically,  this  boundary  is  an  immensely  complex  struc¬ 
ture.  Even  in  the  case  of  the  much  simpler  semiconductor-vacuum  interface 


31-33 


electronic  surface  structure  calculations  have  yet  to  yield  good  quantitative 


35 

agreement  with  microscopic  measurements  .  Undoubtedly,  the  connection  of 
the  microscopic  calculations  with  the  essentially  macroscopic  measurements 
of  capacitance  and  conductance,  which  are  frequently  made  in  studying 
semiconductor-insulator  interfaces  using  the  MOS  structure,  is  even  further 
away.  Although  these  measurements  are  usually  expressed  in  terms  of  inter¬ 
face  state  densities,  it  is  clear  that  there  is  actually  no  real  correspond¬ 
ence  between  such  measurements  and  microscopic  quantities.  We  believe  that 
the  approach  presented  in  this  paper,  which  depends  crucially  on  the  new 
(missing)  semiconduction  boundary  conditions,  affords  a  purely  macroscopic 
quantitative  description  of  semiconductor  interfaces  that  can  be  used  to 
advantage  in  categorizing  semiconductor  interfaces.  More  specifically,  as 
already  noted,  in  this  approach  the  semiconductor  interface  is  characterized 
by  a  set  of  macroscopic  material  surface  coefficients  which  are  to  be 
measured.  In  the  case  of  the  semiconductor- insulator  junction  MOS  quasi¬ 
static  C-V  measurements  provide  a  means  of  evaluating  these  coefficients,  as 
is  exhibited  in  the  last  section  for  the  Si-Si02  interface.  Among  other 
things,  a  knowledge  of  the  values  of  these  material  surface  coefficients  can 
provide  a  detailed  quantitative  understanding  of  the  influence  of  the  surface 
on  the  performance  characteristics  of  devices.  For  example,  in  the  case  of 
a  MOSFET  the  range  in  which  the  values  of  the  surface  coefficients  must  remain 

with  aging  in  order  to  satisfy  some  circuit  performance  criterion  can  be 
36 

established 

2.  Hie  Macroscopic  Model  of  the  Semiconductor 

As  indicated  in  the  Introduction  the  macroscopic  model  of  the  semicon¬ 
ductor  consists  of  three  well-defined  interacting  continua,  which  are  present 
at  each  point  of  space  and  time.  The  three  continua  are  defined  as  follows: 


J 


1.  The  (combined)  lattice  continuum,  denoted  by  the  superscript  i, 

may  have  either  a  positive  or  negative  charge  density  corre- 

37 

sponding  to  the  ionized  impurity  density  .  This  continuum 
does  not  move  and  contains  the  macroscopic  elements  that 
account  for  polarization38. 

2.  The  conduction  electronic  continuum,  denoted  by  the  super¬ 
script  e,  corresponds  to  the  conduction  band  electrons  and  is, 

of  course,  negatively  charged.  This  continuum  is  an  inertialess, 
conducting,  compressible  fluid  that  experiences  a  force  of  re¬ 
sistance  from  its  motion  with  respect  to  the  lattice  continuum. 

3.  The  hole  continuum,  denoted  by  the  superscript  h,  corresponds 
to  the  absence  of  valence  band  electrons  and  is  positively 
charged.  This  continuum  also  is  an  inertialess,  conducting, 
compressible  fluid  that  experiences  a  force  of  resistance  from 
its  motion  with  respect  to  the  lattice  continuum. 

The  Cartesian  components  of  points  in  the  lattice  continuum  are  denoted 
by  x^  (i = 1,2,3),  which,  of  course,  denotes  the  components  of  the  same 
point  in  the  conduction  electronic  and  hole  fluids  also.  The  charge  den¬ 
sities39  associated  with  the  three  continua  are  denoted  p^,  pe  and  p*1, 
respectively,  and  are  functions  of  space  x^  and  time  t  in  general.  The 

three  continua  are  permitted  to  exchange  charge  with  one  another  through 

i  ©  h 

defined  charge  source  densities  labeled  y  ,  y  and  y  .  These  allow  for  the 
possibilities  of  bulk  generation/ recombination  of  carriers  by  means  of 
either  changes  in  the  degree  of  ionization  of  impurities  or  electron-hole 
pair  generation  or  recombination.  Ihen,  in  order  to  satisfy  the  conserva¬ 
tion  of  total  macroscopic  charge  at  each  macroscopic  point,  we  must  have 


(2.1) 


7. 


Now,  each  conducting  fluid  interacts  with  the  lattice  continuum  at  the 

macroscopic  point  x^  by  means  of  defined  material  electric  fields  designated 

6  ll  6  6 

E  and  E  ,  respectively.  These  cause  equal  and  opposite  forces  ±  p  E  and 

±pV  to  be  exerted  between  the  lattice  continuum  and  each  of  the  respective 

'X' 

6  6  h  h 

conducting  fluids.  The  macroscopic  forces  -p  E  and  -p  E  exerted  on  the 

lattice  continuum  by  the  fluids  are  assumed  to  be  sufficiently  small  that 

40 

any  motion  they  tend  to  produce  may  be  ignored  for  our  purposes  .  Each 

conducting  fluid  interacts  with  neighboring  elements  of  the  same  fluid  by 

e  h 

means  of  pressure  forces  labeled  p  and  p  ,  respectively,  which  act  on  the 
surfaces  of  separation  between  elements  of  the  respective  fluids.  The  quasi¬ 
static  Maxwell  electric  field  E  exerts  forces  on  all  elements  of  charge,  but  we  note 
that  since  the  lattice  continuum  does  not  move, only  those  forces  exerted  on  the 
conduction  electronic  and  hole  fluids  are  of  any  importance  in  this  work. 


3.  The  Macroscopic  Equations  for  the  Semiconductor 

Since  the  two  conducting  fluids  have  different  velocities  yet  occupy 
the  same  region  of  space  as  the  lattice  continuum  which  does  not  move,  it  is 
advantageous  to  write  the  integral  forms  of  the  balance  equations  with  respect 
to  a  stationary  element  of  volume.  Accordingly,  we  write  the  charge  balance 
equations  for  the  lattice  (impurity),  the  conduction-electronic  and  the  hole 
continua  in  the  respective  integral  forms 


at  J  pldv  =  J  yldv’ 

V  V 

at  I  pSdv  +  J  5  ’  ^epe<as  =  I  ^dv' 

VS  V 

f  phdV  +  [  n  •  vhphdS  =  f  yhdV  , 
VS  V 


(3.1) 


(3.2) 


(3.3) 


where  n  is  the  outwardly  directed  unit  normal  to  the  surface  S  enclosing  the 

6  ll 

arbitrary  volume  V  fixed  in  space.  The  vectors  v  and  v  denote  the  veloci¬ 
ties  of  the  conduction  electronic  and  hole  continua,  respectively.  The  local 
differential  forms  resulting  from  (3.1)  -  (3.3)  are  obtained  in  the  usual  way 
by  employing  the  divergence  theorem  to  convert  the  surface  integrals  to 
volume  integrals  and  by  employing  the  fact  that  the  volume  V  is  fixed  in 
space  and  is  arbitrary,  with  the  result 


dp  Vat  =  y1 , 

(3.4) 

ape/3t  +  V  *  (peve)  =  ye  , 

(3.5) 

dph/dt  +  V  *  (p hvh)  =  yh  . 

(3.6) 

Clearly,  the  total  charge,  p,  and  the  actual  current  density,  J,  at  any 

point  are  given  by 

e  h  i 

P  =  P  +  P  +  P  > 

(3.7) 

,  e  e  .  h  h  e  ,  h 

J  =  p  v  +  p  v  =  J  +  J  , 

(3.8) 

where  J6  and  are  the  conduction-electronic  and 

<v  /v 

hole  current  densities, 

respectively.  Adding  (3.1)  -  (3.3)  and  employing 

(2.1),  (3.7)  and  (3.8), 

we  obtain 

V  s 

(3.9) 

which  is  the  usual  integral  form  of  the  conservation  of  total  electric  charge. 
From  (3.9),  with  the  aid  of  the  divergence  theorem  and  the  arbitrariness  of  V, 
we  obtain 

9p/f)t  +  V  •  J=  0  , 


(3.10) 


9. 


which  is  the  usual  local  differential  form  of  the  total  charge  balance  equa¬ 
tion.  Note  that  Eq. (3.10)  may  equally  readily  have  been  obtained  by  adding 
(3.4)  -  (3.6)  and  employing  (2.1),  (3.7)  and  (3.8). 

In  this  work  we  need  the  usual  charge  and  circulation  equations  of 
electrostatics,  which  take  the  respective  integral  forms 

J  n  •  DdS  =  J  p  dV  ,  (3.11) 

S~  ~  V 

<f  E  •  dr  =  0  ,  (3.12) 

J  ~  7 

C 

where  D  is  the  electric  displacement  vector,  E  is  the  Maxwell  electric  field 
vector,  S  is  an  arbitrary  surface  enclosing  the  volume  V  and  C  is  an  arbitrary 
closed  curve.  From  (3.11),  with  the  aid  of  the  divergence  theorem  and  the 
arbitrariness  of  V,  we  obtain 

V  *  D  =  p  ,  (3. 13) 

which  is  the  local  differential  form  of  the  charge  equation  of  electrostatics. 
Furthermore,  Eq. (3.12),  with  the  aid  of  Stokes  theorem  and  the  arbitrariness 
of  the  area  enclosed  by  C,  enables  us  to  define  the  usual  scalar  electric 
potential  cp  such  that 

E  =  -  Vcp  .  (3.14) 

We  now  note  that  since  the  lattice  continuum  does  not  move  the  balances 
of  linear  momentum  and  mass  of  the  lattice  continuum  may  be  ignored  in  this 
work  in  which  the  lattice  continuum  simply  serves  as  the  entity  that  provides 
the  reactions  to  the  actions  of  the  two  conducting  fluids.  Furthermore, 
since  the  conducting  fluids  have  been  defined  as  being  massless,  the  balance 
of  mass  of  these  constituents  is  not  needed.  The  equations  of  the  balance 
of  linear  momentum  for  the  conduction-electronic  and  hole  continua  take  the 


respective  integral  forms 


10 


-  J  npedS  +  J  pe(E+Ee)  dv  =  0  ,  (3.15) 

S  V 

-J  nphdS  +  J  ph(E  +  Eh)  dV=  0  .  (3.16) 

S  V 

From  (3.15)  and  (3.16),  with  the  aid  of  the  divergence  theorem,  the  arbitrari¬ 
ness  of  V  and  (3.14),  we  obtain  the  local  differential  forms  of  the  semicon¬ 
duction  field  equations  for  the  respective  continua 

"  XP6  "  +  =  0  j  0.17) 

-  Vph  -  phVcp  +  pV*  =  0  .  (3.18) 

At  this  point  we  note  that  the  governing  differential  equations  consist¬ 
ing  of  (3.10),  (3.13),  any  two  of  (3.4)  -  (3.6),  (3.17)  and  (3.18),  with 
(2.1),  (3.7),  (3.8)  and  (3.14),  constitute  an  underdetermined  system,  i.e., 
there  are  more  dependent  variables  [31,  (10  scalars  and  7  vectors)]  than  equa¬ 
tions  (18) .  As  usual  in  any  macroscopic  description  of  this  nature,  materi¬ 
ally  descriptive  constitutive  equations  are  required  in  order  to  obtain  a 
determinate  system.  Constitutive  equations  consistent  with  the  model  are 
obtained  by  writing  the  equation  of  the  conservation  of  energy  for  the 
combined  material  continuum,  from  which  the  first  law  of  thermodynamics  is 
obtained,  and  employing  the  second  law  of  thermodynamics.  However,  in  order 
not  to  obscure  the  main  purpose  of  this  work,  which  is  to  indicate  the 
importance  of  the  boundary  conditions  as  well  as  the  differential  equations 
in  the  macroscopic  description  of  semiconductors,  we  give  the  derivation  of 
consistent  constitutive  equations  from  the  thermodynamics  in  Appendix  A 
and  simply  present  the  resulting  constitutive  relations  here.  In  accordance 
with  Appendix  A  the  constitutive  equations  must  take  the  following  functional 


forms: 
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pS  =  pe(pe,T),  ph  =  ph(ph,T), 

(3.19) 

£e  =  Ee<pV>S,T>>  Eh  =  Eh(Ph,vh,E,T), 

(3.20) 

5e  =  Se<E,T)  , 

(3.21) 

e  e ,  e  h  e  h  „ 

y  =y  (p  ,p  ,  v  ,  v  ,E,T)  , 

h  h.  e  h  e  h  _ 

Y  -Y  <P  ,P  ,v  ,E,T)  , 

(3.22) 

where  T  is  the  (here  uniform)  absolute  temperature.  As  in  all  macroscopic 
theories,  e.g.,  electromagnetism  and  elasticity,  the  functional  forms  and 
the  values  of  the  material  coefficients  (dielectric  constants,  elastic 
constants,  etc.)  are  to  be  determined  from  experiment.  Any  functional  form 
depending  on  the  variables  shown  in  (3.19)  -  (3.22)  will  be  consistent  with 
the  model  provided  that  (A7)  is  satisfied.  For  example,  the  important  forms 
usually  chosen  for  (3.20)  are 


_e  e  .  e  h  h.  h 
E  =  v  /)i  ,  E  =-v/n  , 


(3.23) 


where  p,6  and  p,*1  are  the  measured  mobilities  of  the  conduction-electronic  and 

hole  fluids,  respectively.  At  high  electric  fields,  it  is  well  known  that 
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E  dependent  terms  become  necessary  .  Similarly,  the  form  usually  chosen 
for  (3.21)  is 


D 


e 


E 

s~ 


> 


(3.24) 


where  e  is  the  permittivity  for  the  semiconductor, 
s 

Before  proceeding,  we  note  that  a  significant  simplification  in  (3.17) 

and  (3.18)  may  be  achieved  by  employing  certain  relations  between  the 

0  h 

gradients  of  the  pressures  p  and  p  and  the  gradients  of  the  well-defined 

6  6  0  }j  h  h  « 

chemical  potentials  cp  =  cp  (p  ,T)  and  cp  *cp  (p  ,T),  which  are  derived  in 

Appendix  A  and  are  of  the  form 
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(l/p®)7p®  -  7cp®  ,  (l/ph)7ph  -  7cph  .  (3.25) 

The  substitution  of  (3.25)  into  (3.17)  and  (3.18)  yields 

-7(cp®  +  cp)  +  E®=0,  (3.26) 

-7(cph  +  cp)  +  Eh  =  0  ,  (3.27) 


which,  as  we  shall  see,  are  very  useful  forms  of  the  differential  equations 
for  the  semiconducting  fluids.  Substituting  from  (3.23)  for  the  respective 
fluids  into  (3.26)  and  (3.27),  respectively,  and  employing  (3.14),  the 
definitions  in  (3.8),  the  fact  that  cp®  =  cp®  (p®,  T)  and  cp*1  =  cp*1  (p*1, T)  and  the 
assumption  of  uniform  temperature,  we  obtain 

J®  =-p.epeE-De7pe  ,  (3.28) 

3*  =  (A>V°V\  (3-29> 

where 


_e  e  e,  e  e 
D  =  -  p,  p  Bcp  /dp 


h  .  h  h  h  h 
D  =  +  p,  p  dcp  /dp  , 


(3.30) 


.  43 

are  the  diffusion  coefficients.  Equations  (3.28)  and  (3.29)  are  the  standard 

diffusion-drift  current  equations  for  the  conduction-electronic  and  hole 

44 

fluids,  respectively,  and  (3.30)  contains  familiar  expressions  also  . 
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Appropriate  forms  for  (3.22)  are  well  known  . 

At  this  point  we  note  that  we  now  have  a  determinate  theory,  which  by 
appropriate  substitution  can  readily  be  reduced  to  10  equations  in  the  10  de¬ 
pendent  variables,  cp,  p,  p®,  p*1,  Jand  J®.  The  10  equations  are  any  two  of 

(3.4)-(3.6),  (3.10),  (3.13)  and  the  three  each  of  either  (3.26)  and  (3.27)  or 
46  47 

(3.28)  and  (3.29)  *  .  In  order  to  have  a  complete  field  theory,  the 

boundary  (or  jump)  conditions  across  surfaces  have  to  be  adjoined  to  the  sys¬ 
tem  of  differential  equations.  However,  before  we  discuss  the  boundary 


13 


conditions  we  note  that  the  particularly  useful  differential  forms  for  the 
semiconducting  fluids  are  equivalent  to  the  integral  forms 

-J  n(cpe+cp)dS  +  J  E®dV«0,  (3.31) 

-J  n(cph  +  cp)dS  +  J  EhdV-0,  (3.32) 

because  the  local  differential  forms  in  (3.26)  and  (3.27)  can  readily  be 
obtained  from  the  respective  integral  forms  in  (3.31)  and  (3.32)  when  the 
field  variables  are  differentiable.  Moreover,  we  now  take  (3.31)  and  (3.32) 
to  be  valid  even  when  the  field  variables  are  not  differentiable  and  (3.26) 
and  (3.27)  cannot  be  obtained,  such  as  across  surfaces  of  discontinuity. 

Of  course,  all  the  other  integral  forms  are  taken  to  hold  across  surfaces 
of  discontinuity  also.  For  purposes  of  obtaining  boundary  conditions  in  this 
work,  the  integral  forms  (3.31)  and  (3.32)  replace  (3.15)  and  (3.16),  respect¬ 
ively,  because  they  are  considerably  simpler  forms  that  result  in  more  con¬ 
venient  boundary  conditions. 

The  pertinent  boundary  (or  jump)  conditions  are  obtained  by  applying 

the  integral  forms  (3.9),  (3.11),  (3.31),  (3.32)  and  any  two  of  (3.1)  -  (3.3) 

to  an  arbitrary  pill  box  region  encompassing  a  portion  of  the  surface  of 

discontinuity  and  taking  the  limit  of  the  region  as  the  volume  shrinks  to 

48 

zero  faster  than  the  area  in  the  usual  way  ,  while  assuming  that  the 
volumetric  densities  either  remain  bounded  or  become  unbounded  in  a  specified 
way.  The  boundary  (or  jump)  condition  obtained  from  (3.12),  which  contains 
a  circulation  integral,  is  determined  instead  by  taking  the  circulation  around 
an  arbitrary  area  normal  to  and  intersecting  a  portion  of  the  surface  of  dis¬ 
continuity  and  taking  the  limit  as  the  area  enclosed  by  the  circuit  collapses 

48 

to  the  surface  of  discontinuity  in  the  usual  way  .  The  jump  conditions 
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obtained  from  any  of  (3.1)  -  (3.3)  are  not  of  interest  in  the  present  work 
(because  only  total  surface  charge  is  of  interest)  and  consequently,  are  not 
obtained  here.  From  (3.11)  in  the  above  mentioned  way  we  obtain  the  well- 
known  condition 

n  •  [D]  =  ct  ,  (3.33) 

where  we  have  introduced  the  conventional  notation  [CJ  for  C+  -  C  ,  n  denotes 
the  unit  normal  directed  from  the  -  to  the  +  side  of  the  surface  of  discon¬ 
tinuity  and  a  is  the  surface  charge  density,  which  is  defined  by 

f  a  dS  =  lim  j  p  dV  .  (3.34) 

i  V-0  i 


Similarly,  from  (3.9)  with  the  aid  of  the  time  derivatives  of  (3.33)  and  (3.34), 
we  obtain 

n  •  [dD/dt  +  J]  *  0  .  (3.35 

From  the  circulation  integral  (3.12)  in  the  above  mentioned  usual  way,  we 
obtain  the  well-known  condition  of  electrostatics 


n  X  [E]  =  0  , 


(3.36) 
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which,  with  the  aid  of  (3.14),  enables  us  to  obtain  the  boundary  condition 

_[C<5]  “  0  •  (3.37) 

The  application  of  (3.31)  to  the  arbitrary  limiting  pill  box  region 
in  the  aforementioned  manner  yields 


where 


6  G 

n[c p  +  cp]  ■  f  n  , 


[  fen  dS  -  lim  f  Ee  dV  , 
«l  ~  V-»0  •>  ~ 


(3.38) 


(3.39) 


and  we  note  that  the  force  per  unit  charge  E  exerted  by  the  lattice  continuum 


on  the  conduction  electronic  fluid  at  the  interface  must  be  permitted  to  become 


unbounded  in  the  manner  indicated  in  (3.39)  because  there  must  be  a  force 
exerted  by  the  lattice  continuum  on  the  conduction  electronic  fluid  at 
the  boundary  which  prevents  the  electrons  from  leaving  the  semiconductor. 
Substituting  from  (3.37)  into  (3.38)  and  omitting  the  n,  which  is  unneces¬ 
sary,  we  obtain 

(3.40) 

which  states  that  the  jump  in  the  chemical  potential  of  the  conduction  electronic 
fluid  across  a  surface  of  discontinuity  is  equal  to  the  restraining  surface 
interaction  force  exerted  by  the  lattice.  Clearly,  from  (3. 20)^ and  (3.39), 
the  constitutive  equation  for  f  must  be  of  the  form 

fe=fe(pe,ve,E,T)  .  (3.41) 

In  a  similar  way,  from  (3.32)  and  (3.37),  we  obtain 

[cp*1]  “  fk  ,  (3.42) 

where 

f  fhndS=  lim  f  EedV,  (3.43) 

<1  ~  Q  J  ~ 

S  V 

and  from  (3.20)2and  (3.43)  we  must  have 

fh.£h(ph,vh,E,I).  (3.44) 

Equations  (3.40)  and  (3.42),  with  (3.41),  (3.44)  and  the  functional  depend- 
0  ll  0  1^ 

ences  of  cp  and  cp  ,  respectively,  on  p  and  p  are  the  previously  discussed 
required  semiconduction  boundary  conditions,  the  existence  of  which  along 
with  the  application  in  the  next  section  constitute  the  main  point  of  this 
work.  Appropriate  specific  forms  for  (3.41) and  (3.44)  are  discussed  in  the 
application  in  the  next  section  and  expressions  for 

cp6  *  cp6  (pe,T) ,  cp*1  ■  cp*1  (p**,T)  (3.45) 


are  discussed  in  Appendix  B. 
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4.  Application  to  Si-Sio2  Interfaces 

As  noted  in  the  Introduction,  the  purely  macroscopic  equations  for  the 
semiconductor,  which  have  been  developed  here,  are  particularly  well-suited 
for  the  description  of  the  influence  of  the  properties  of  the  surface  of 
the  semiconductor  on  gross  sample  behavior.  The  MOS  quasi-static  C-V  experi¬ 
ment  provides  a  cohvenient  means  of  demonstrating  some  of  the  advantages  of  the 
purely  macroscopic  description.  As  already  mentioned,  this  experiment 
provides  a  means  of  evaluating  some  of  the  macroscopic  material  surface 
coefficients  that  occur  in  expressions  to  be  written  for  fe  and  f*1. 

The  values  of  these  coefficients  for  a  particular  semiconductor 
surface  are  obtained  by  solving  a  straightforward  one-dimensional  problem 
corresponding  to  the  experimental  setup  and  deriving  an  expression  for  an 
observable  consequence  in  terms  of  the  unknown  material  surface  coefficients, 
from  which  the  values  of  the  coefficients  are  determined  by  comparison  with 
measurements.  In  order  that  our  model  of  the  actual  MOS  quasi-static  C-V 
experiment  be  valid,  the  following  requirements  must  be  satisfied^: 

1.  The  geometry  should  be  such  that  the  assumption  of  one¬ 
dimensional  dependence  holds^. 

2.  The  thickness  of  the  semiconductor  should  be  sufficiently  large 
that,  one  end  can  be  assumed  (in  the  problem)  to  be  at  infinity. 


3.  There  should  be  no  fixed  or  mobile  oxide  charges. 
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4.  The  static  form  of  the  equations  should  be  applicable' 

5.  The  semiconductor  properties  should  be  macroscopically  uniform. 

6.  The  oxide  capacitance  should  not  be  much  smaller  than  other 

capacitances  in  order  that  the  surface  behavior  be  detectable 
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by  a  total  capacitance  measurement 


A  schematic  of  the  modelled  MOS  structure  with  associated  coordinate  system 
is  shown  in  Pig.l. 


Under  the  assumed  static  conditions  we  have 


eh.  e  h  i 
v  -v  -0,  y  =  y  -V  -0, 


(4.1) 


and  the  dynamic  equations  consisting  of  (3.5),  (3.6)  and  (3.10)  are  satis¬ 
fied  identically.  Hie  nontrivial  equations  remaining  are  (3.13),  (3.28) 
and  (3.29)  in  the  absence  of  Je  and  J*1.  Since  in  the  problem  under  consid- 
eration  we  have  one-dimensional,  i.e.,  x-,  dependence  only,  from  (3.13), 
with  (3.7),  (3.14)  and  (3.24),  we  have 


2  .  2  .  .  .  e  h  i.  . 

d  cp/dx  =  -  (1/e  )  (p  +p  +  p  )  ,  O^x, 


and  (3.28)  and  (3.29)  take  the  respective  forms 


pBdcp/dx=  (De/p,e)dpS/dx  ,  0  <.  x, 


(4.2) 


p^dcp/dx  =  -  (Dh/y,h)  dph/dx  ,  0  <.  x  . 


Since  the  charge  density  vanishes  in  the  oxide,  we  have 


(4.3) 

(4.4) 


d^cp/dx^  =  0  ,  SxsO  .  (4.5) 

As  a  consequence  of  condition  2  above,  we  have 

p  -♦  0  ,  cp-*0,  dcp/dx-*0,  as  x-»oo.  (4,6) 

The  boundary  condition  to  be  satisfied  at  the  interface  between  the  oxide 
and  the  metal  is 


'p  =  'Po"'Pmos>  (4-7) 

where  <p„__  is  due  to  the  potential  jumps  arising  from  possible  double  layers 

at  the  metal-oxide  interface  and  the  metal-semiconductor  interface  and 

a>  is  the  applied  (gate)  voltage  relative  to  a  metal  "ground"  plate. 

~o 

From  (3.37),  (3.40)  and  (3.42)  we  see  that  the  additional  boundary  conditions 


to  be  satisfied  at  the  interface  between  the  oxide  and  the  semiconductor  are 


18. 


j<pe]  =  fe  (pe,  dcp/dx]  ox  ,  dcp/dx I  sem  ,  T)  ,  at  X  *  0  , 

[cph]  =  f*  (ph,  dcp/dx|ox  ,  dcp/dx|  gem  ,  T)  ,  at  x  =  0  , 

where  the  electric  field  terms  are  evaluated  on  either  side  of  and 
adjacent  to  the  oxide-semiconductor  interface,  and  we  note  that 


(4.9) 

(4.10) 

immediately 


*o*(0)  =  ®se»10’  *'V  (4-11» 
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where  cp  is  the  surface  potential  of  the  semiconductor  .  Integrating  (4.3) 
s 

and  (4.4)  and  employing  the  boundary  conditions  (4.6),  we  obtain 

e  h  . 

cp  =  —  iog  pe/p®=  -  log  p"/p*  ,  0  ^  x,  (4.12) 

0  h 

where  p  and  p  are  the  constant  charge  densities  of  the  respective  fluids 

CO  CD 

at  infinity.  Prom  (4.6)^,  condition  5  and  (3.7)  we  have 

P®  +  P®+Pi  =  0-  (4'13) 


0  ll 

Solving  (4.12)  for  p  and  p  in  terms  of  cp,  substituting  in  (4.2)  and 
employing  (4.13),  we  obtain  a  single  nonlinear  differential  equation  in  the 
single  dependent  variable  cp,  a  first  integral  of  which  may  be  written  in 
the  form 


,  r  _e  ,  e  e  e  .  ,  h  .  h  .  h  h  ._l/2 


(4.14) 


where  the  integration  constant  has  been  found  from  the  boundary  condition 
ait  infinity  given  in  (4.6)^.  The  sign  of  dco/dx  is  undetermined  at  this 
stage.  Equation  (4.14)  may  be  reduced  to  a  quadrature  immediately.  In 
order  to  evaluate  the  integration  constant  and  fix  the  presently  unknown 

j 

sign  in  (4.14),  we  must  esploy  the  heretofore  unknown  (missing)  boundary 
conditions  in  (4.9)  and  (4.10).  Since  these  boundary  conditions  involve 


J 
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the  chemical  potentials  cp  and  cp  ,  we  digress  with  a  derivation  of  some  simple 
e  li 

expressions  for  cp  and  cp  in  the  static  case.  By  virtue  of  (3.30),  (4.3) 
and  (4.4)  can  be  written  in  the  form 

dcp/dx  +  dcp6/dx  =  0,  OSx  ,  (4.15) 

dcp/dx  +  dcp  /dx  =  0,  0  £  x  ,  (4. 16) 

which  are  just  the  one-dimensional  versions  of  (3.26)  and  (3.27)  in  the 
e  h 

absence  of  E  and  E  ,  which  vanish  when  the  velocities  of  the  conduction 
electronic  and  hole  fluids  vanish. 

Clearly,  Eqs. (4.15)  and  (4.16),  which  are  the  same  as  (4.3)  and  (4.4), 
can  be  integrated  immediately  to  obtain 

cp+cp6  =  cp®,  Osx,  (4.17) 

(p+<ph*cp^,  OSx,  (4.18) 


0  h 

where  cp  and  cp  are  the  constant  values  of  the  chemical  potentials  away  from 

1  CO  1 00 
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the  space  charge  region,  and  we  have  employed  (4.6)  .  Now,  the  left-hand 

sides  of  (4.9)  and  (4.10)  are 

h 


r  e,  e  e  .  h,  h 
[cp  J  =  cp  -cp  ,  [cp  J  =  cp  -cp  , 

~  i  sem  ^ox  ’  ~  ~  rsem  ox  ’ 


(4.19) 


@  h 

into  which  we  can  substitute  from  (4.17)  and  (4.18)  for  cp  and  cp 

sem  sem 

])  e 

Subtracting  [cp  ]  from  [cp  3 ,  we  obtain, 


,  e,  h  e  h  e  .  .  h, 

i<P  i  =  9ox-9ox-cpa5  +  cpa3+ 19  L> 


(4.20) 


eh 

where,  of  course,  cp  and  cp  are  the  constant  chemical  potentials  in  the 

OX  ox 

oxide.  The  four  constants  appearing  in  (4.20)  are  presumed  known  and, 

e  h 

therefore,  in  the  static  case  [cp  3  and  [cp  3  and,  of  course,  from  (4,9)  and 


(4.10)  f  and  f  are  not  independent  quantities.  Consequently,  in  the 


static  case  (4.9)  and  (4.10)  are  not  independent  boundary  conditions,  which 
is  not  surprising  since  only  one  integration  constant  remains  undetermined. 


Since  fe  and  r  differ  only  by  a  constant  and  the  functional  dependences  for 
G  ll  ^ 

f  and  f  are  of  the  forms  shown  in  (4.9)  and  (4.10),  it  is  clear  from  (4.12) 
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that  the  functional  dependences  of  f  and  r  on  p  and  p  must  be  through 

Q  0  ll  ll 

the  respective  variables  log  p  /p  and  log  p  /p  because  the  conditions  (4.9), 

0  h 

(4.10)  and  (4.20)  hold  for  arbitrary  p  and  p  .  In  view  of  the  foregoing, 
from  (4.9),  (4.11)^  (4.17)  and  (4.19),  the  one  required  boundary  condition 
at  the  semiconductor-oxide  interface  in  the  static  case  must  be  of  the  form 

cpG-cpe  -cp  =fe(cp  ,dcp/dx|  ,  dcp/dxj  ,T)  .  (4.21) 

T°°  Tox  Ys  Ts’  1  ox  *  r  '  sem' 

Equation  (4.21)  provides  a  relation  between  the  surface  potential  cp  and  the 
~  s 

electric  fields  on  either  side  of  the  interface  and,  as  noted  earlier, 

contains  material  surface  coefficients,  which  are  to  be  measured  and  are 

characteristic  of  the  particular  interface.  When  the  material  surface 

coefficients  are  known,  Eq. (4.21)  provides  the  one  independent  condition 

which  enables  the  one  remaining  integration  constant  to  be  found  and,  thus, 

the  solution  to  the  posed  problem  to  be  obtained. 

In  the  case  of  the  static  one-dimensional  problem  being  treated  here 

calculations  are  greatly  simplified  owing  to  (4.14),  which  may  be  evaluated 

at  x  =  0  to  obtain  dcp/dxl  explicitly  as  a  function  of  cp  .  By  satisfying 

1  sem  T  s 

(4.5),  (4.7)  and  (4.11)  for  the  oxide  portion  of  the  problem,  we  find 

dcP/dx)ox=  %-%  +  Vms)/z>  (4’22) 

and,  hence,  (4.21),  with  (4.14)  and  (4.22)  provides  a  direct  relation  between 
cpQ  and  cpg  at  constant  temperature.  The  aforementioned  simplification  in 
the  calculations  also  causes  a  weakness  in  our  ability  to  obtain  definitive 
values  of  the  material  surface  response  coefficients  from  MOS  quasi-static 
one -dimensional  C-V  measurements  and  to  distinguish  conclusively  between 
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field  and  surface  potential  (charge  density)  terms  in  the  expression 
0 

for  f  .  In  this  paper,  where  we  intend  only  to  illustrate  the  completely 
macroscopic  approach  to  semiconductor  interfaces,  we  are  content  merely 
to  provide  reasonably  convincing  (but  not  definitive)  evidence  for  the 
forms  chosen  for  (4.21). 

0 

In  establishing  a  form  for  f  and  determining  the  values  of  the 

attendant  material  surface  coefficients  that  arise,  experimental  data  is 

utilized.  In  the  quasi-static  C-V  experiment  the  total  capacitance  is 

measured  directly  as  a  function  of  the  applied  voltage.  The  Berglund  inte- 
.  56 

gration  of  the  capacitance  may  be  used  to  give  reasonably  accurate  values 
,  57 

for  cps  vs  cp^  over  certain  ranges  of  cpg  .  The  aforementioned  simplification 

of  (4.21)  resulting  from  the  existence  of  (4. 14)  makes  this  (m  vs  qj  )  data 

Ts  Yo 
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preferable  .  It  is  also  to  be  recommended  by  the  fact  that  this  data  is 
a  direct  measure  of  semiconductor  surface  response,  whereas  total  capacitance 
is  not.  At  this  point  it  should  perhaps  be  mentioned  that  the  procedure  of 
using  the  C-V  data  to  generate  [- (l/q)da/dcpsl ,  which  is  commonly  referred 
to  as  N  ,  the  surface  state  density,  is  unnecessary  insofar  as  the  determin- 
ation  of  macroscopic  surface  coefficients  is  concerned.  However,  in  a  number 
of  cases  we  have  used  this  data  since  it  is  the  form  usually  presented  in  the 
literature. 

.  .  G 

Now,  an  explicit  representation  of  f  in  (4.21)  is  found  by  choosing  a 

suitable  simple  expression  in  the  variables  shown  and  obtaining  the  material 

coefficients  in  the  expression  from  experimental  c  vs  cp^  or  dc/dcps  vs  cpg  data. 

59 

Hie  wide  range  of  variation  in  measured  surface  response  seems  to  indicate 

that  the  possibility  of  obtaining  one  form  for  fe  for  all  Si-SiC^  interfaces 

is  unrealistic.  However,  for  annealed  MOS  structures,  interface  response 

59  60 

is  found  to  be  qualitatively  similar  and  hence  for  these  interfaces  a  single 
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form  for  f  appears  to  be  appropriate.  The  expression  we  choose  for  f  is 


+  A  (’ T.m  )  +  A_  (T>C0  +  A.(T)cd^+.. 


f  =  A1(T'V  dir  +A2(T’CPS)  d£  +  A3  (T>cps  +  A4  (T)cp‘ 

ox  sem 


(4.23) 


where  all  coefficients  are  functions  of  temperature,  which  is  constant  in  the 
present  work.  We  note  that  the  electric  field  terms  must  be  present  in  (4.23) 
because  the  dependence  of  c pg  on  cpo  arises  only  through  the  dcp/dx | ^  term  and 
unless  a  dependence  on  cpQ  appears  in  (4.23),  (4.21)  reduces  to  an  equation  for 
cps  independent  of  the  external  conditions.  Only  terms  linear  in  dcp/dx  on  each 
side  of  the  interface  are  retained  in  (4.23)  because  the  electric  field  is 
the  force  per  unit  charge.  Further  specification  of  (4.23),  as  to  which  terms 
are  to  be  kept  and  the  values  of  the  coefficients  A^,  is  accomplished  by  use 
of  a  regression  analysis  (including  plots  of  residuals)  of  the  experimental 

data.  The  goal  here  is  to  achieve  reasonable  accuracy  in  the  representation 
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of  f  while  employing  the  minimum  number  of  terms  .  However,  it  is  important 

to  realize  that,  ultimately,  the  level  of  accuracy  required  in  the  representa- 
.  6 

txon  of  f  is  dictated  by  the  particular  problem  under  consideration.  It  is 

found  that  although  the  surface  coefficients  A,  and  A_  can  depend  on  cp  ,  no 

1  &  Ts 

dependence  on  cpg  is  required  in  order  to  obtain  reasonable  agreement  with  the 
experimental  data  considered  in  this  work.  However,  such  a  dependence  can 
improve  the  accuracy  of  the  representation.  The  remaining  terms  in  the  expres- 

.  G 

sxon  for  f  are  a  power  series  in  cpg.  The  number  of  terms  retained  in  any 
polynomial  approximation  depends  on  the  accuracy  desired.  We  have  found  that 
typically  no  more  than  three  terms  in  the  power  series  are  necessary  to  give 
adequate  agreement. 

For  the  purpose  of  obtaining  the  material  surface  coefficients  Ai  from 
experimental  data  of  a  vs  cpg  or  da/dcpg  vs  cpg,  the  solution  can  be  put  in  a 
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particularly  useful  form  by  employing  the  constitutive  relation  for  the 
oxide 


D  =  -  e  dcp/dx  , 
ox  ox  Y  * 

where  eQX  is  the  permeability  of  the  oxide.  By  substituting  from 
and  (4.24)  into  (3.33)  we  obtain  an  expression  for  dcp/dx Jqx  which 
into  (4.23).  The  resulting  expression  for  fe  is  substituted  into 
which  for  A^  /  0  we  obtain 


CT  =  B 

o 


+  B2cps  +  B3cps  + 


> 


(4.24) 
(3.24) 

is  inserted 
(4.21) , from 

(4.25) 


where 
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B2  =  - 
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(1  +  a3) 


i  :>3  , 


(4.26) 


and  we  have  taken  the  liberty  of  omitting  the  temperature  dependences. 

Although  (4.25)  is  a  particularly  convenient  form  for  use  in  this  work,  it 

should  be  noted  that  the  A^  are  to  be  regarded  as  fundamental  here  and  not  the 
6  3 

derived  b”  .  In  existing  work  on  semiconductor  interfaces  the  coefficient  B 
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is  often  called  the  surface  state  charge  .  We  also  remark  that  since  N 

data  is  proportional  to  da/dcp  it  does  not  provide  a  measure  of  B  .  In 

s  o 

Fig.  2  some  experimental  dc/dcp  vs  cp  data^"*  is  presented  along  with  a  reasonably 

s  s 

accurate  fit  to  that  data  using  the  cpg  derivative  of  (4.25).  The  non-zero  B^ 
and  A.  coefficients  determined  from  this  fit  for  this  particular  Si-Si02 
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interface  are  also  shown 


ss 
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An  equation  similar  to  (4.25)  has  been  used  elsewhere  as  the  addi¬ 


tional  boundary  condition  required  for  the  solution  of  a  boundary  value  problem 
to  be  obtained.  However,  as  noted  in  the  Introduction,  assumed  conditions,  such 


24. 


as  that  in  Ref. 23,  are  used  simply  because  something  is  needed  and  are 

not  obtained  from  fundamental  physical  principles  dictated  by  the  field 

6  7 

theory  as  is  done  here  .  Furthermore,  the  assumed  condition  of  Ref. 23 
is  a  static  condition  and  not  dynamic  (as  it  is  used)  as  shown  by  the  deriva¬ 
tion  of  (4.25).  Moreover  the  condition  employed  in  Ref. 23  does  not  contain 
the  term  appearing  in  (4.25).  The  appearance  of  the  term  in  (4.25)  is 
guaranteed  by  the  fact  that  the  terms  linear  in  the  electric  field  must  be 
present  in  (4.23),  as  noted  earlier.  Nevertheless,  it  is  entirely  possible 
that  the  (field)  term  could  be  much  smaller  than  the  other  terns  and, 
hence,  negligible.  It  is  to  be  expected  that  if  the  B^  term  is  ever  import¬ 
ant,  it  would  be  when  very  high  fields  exist  near  the  surface,  which  is  the 
case  when  large  voltages  are  applied  or  when  the  semiconductor  is  heavily 
doped  and  subject  to  moderate  voltages.  This  expectation  is  also  consistent 
with  the  fact  that  high  fields  tend  to  cause  many  quantities  to  be  field 
dependent,  e.g.,  mobilities.  The  regression  studies  appear  to  support 
this  expectation.  However,  as  already  noted,  on  account  of  a  weakness  in 
the  static  one-dimensional  problem  no  definitive  experimental  evidence  of 
the  necessity  of  the  B^  term  can  be  obtained  in  this  way. 

Figure  3  presents  a  comparison  of  fits  of  the  cpg  derivative  of  (4.25),  both 

with  and  without  the  B^  (field)  term,  to  a  particular  data  set®8.  As  was 

typically  the  case  for  the  data  considered  in  this  work,  inclusion  of  the 

field  term  yielded  a  better  fit  with  fewer  terms  needed.  In  addition  in 

strong  accumulation  or  inversion  it  is  the  field  term  which  constitutes  the 

0 

dominant  contribution  to  f  ,  whereas  for  n-type  semiconductors  if  our  sole 
interest  is  in  the  surface  potential  range  -.  55  V  Scpg  <  0. 0  V,  the  field 
term  may  be  accurately  expanded  in  a  few  cpg  polynomial  terms.  We  note  that 
for  the  particular  sample  of  Fig. 3  the  surface  field  dcp/dxl 

sem 

accumulation  is  ~  10®  volts/cm®9. 


in  strong 


The  findings  presented  here  concerning  values  of  material  surface 

coefficients  and  of  the  actual  expression  used  for  f  in  (4.23)  should  be 

regarded  merely  as  illustrative.  Extensive  and  systematic  experimentation 

e 

would  be  needed  to  definitively  establish  a  "best"  form  for  f  for  annealed 

samples.  With  the  aid  of  such  a  form,  a  catalogue  of  the  values  of  surface 

coefficients  along  with  the  sample  preparation  techniques  necessary  to  achieve 

or  alter  those  values  could  be  obtained  and  would  constitute  the  complete 

macroscopic  description  of  the  Si-siC>2  interface.  As  an  illustration. 

Table  I  gives  a  tabulation  of  values  for  material  surface  coefficients  for 

70 

a  specific  Si-SiO_  interface  ,  the  dc/dcp  vs  cp  curves  of  which  are  shown 
6  s  s 

in  Pig. 4,  based  on  the  form  of  f  in  (4.23).  Each  sample  underwent  the  same 
preparation  technique  except  that  different  metals  were  used  for  the  front 
contact  in  each  case. 
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Basic  Semiconductor  Transport  Equations  -  II,”  Solid-State  Electron, 

17,  707  (1974). 

48.  Ref. 9,  Sec. 1.13. 

49.  Although  (3.36)  implies  only  that  [cp]  =  constant,  we  assume  the  constant 
to  be  zero.  This  discounts  the  possible  existence  of  a  double  layer 

at  this  interface.  While  such  things  are  quite  important  for  metal 
surfaces,  we  feel  that  for  semiconductor-insulator  boundaries  such  a 
quantity,  if  not  zero,  is  surely  negligible.  However,  as  with  metal 
boundaries,  it  could  be  included  if  desired. 

50.  All  assumptions  except  condition  6)  may  be  greatly  modified,  however 
since  the  experiment  is  being  used  to  measure  coefficients  defined  by 
the  theory,  the  experiment  can  be  tailored  to  best  meet  the  needs  of 
the  theory,  i.e.,  to  satisfy  conditions  1-5. 

51.  This  requires  that  the  ratios  of  oxide  thickness  and  Debye  length  to 
electrode  width  be  small. 

52.  This  demands  that  the  quasistatic  C-V  experiment  be  done  at  a  low 
enough  frequency  «  10  Hz)  so  that  the  slowest  response  -  that  of  the 
generation  of  minority  carriers  in  the  inversion  layer  -  is  capable  of 
following  the  low  frequency.  See  M.  Kuhn  and  E.H.  Nicollian, 
"Nonequilibrium  Effects  in  Quasistatic  MOS  Measurements,"  J.  Electrochem. 
Soc.,  118.  370  (1970). 

53.  This  requires  that  the  oxide  thickness  not  be  too  great. 

54.  The  assertion  in  (3.37)  makes  (4.11)  unambiguous. 
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55.  In  the  static  problem  no  specification  of  cp  and  cp  is  necessary 
beyond  their  functional  dependences  and  their  values  at  infinity 
and  in  the  oxide.  In  a  dynamic  situation,  where  (4.17)  and  (4.18) 
do  not  hold,  specific  choices  for  cpe  and  cp'1  would  have  to  be  made, 
such  as  those  appearing  in  Appendix  B. 

56.  C.N.  Berglund,  "Surface  States  at  Steam-Grown  Silicon-Silicon  Dioxide 
Interfaces,"  IEEE  Trans.  Electron  Devices,  ED-13,  701  (1966). 
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57.  Me  remark  that  the  low  frequency  C-V  measurement  is  inadequate  for 
probing  semiconductor  surface  response  in  the  extremes  of  surface 
potential,  as  is  well  known.  The  energy  level  ranges  which  correspond 
to  these  surface  potentials  have  been  explored  by  the  method  of  P.V. 

Gray  and  D.M.  Brown,  "Density  of  Si-Si02  Interface  States,"  Appl.  Rhys. 

Lett.,  (3,  31  (1966).  This  method,  however,  is  not  relevant  to  the 
problem  being  treated  in  Sec. 4  because,  within  the  macroscopic  framework 
and  terminology,  it  is  concerned  with  the  temperature  dependence  of  f® 
at  zero  applied  voltage,  whereas  we  are  concerned  with  establishing  the 
functional  dependence  of  f?  at  constant  temperature.  Macroscopically 
speaking,  these  are  distinct  questions. 

58.  Any  quantities  which  can  be  derived  from  cp  vs.  cp  ,  such  as  a  , 

s  To’  sem’ 

a  -(1/q)  da  /dtp  =N  vs.  cp  may,  of  course,  also  be  used.  Another 
met  sem  s  s  s 

measurement  Which  is  perhaps  better  for  the  purposes  of  this  paper  in 
that  it  claims  greater  accuracy  is  the  direct  measurement  of  the  gate 
surface  charge.  See  J.  Koomen,  "The  Measurement  of  Interface  State 
Charge  in  the  MOS  System,"  Solid  State  Electron.,  JL4,  571  (1971). 

59.  B.E.  Deal  and  R.R.  Razouk,  "Structural  and  Electrical  Investigation  of 
the  Si-Si02  Interface  in  Thermally  Oxidized  Silicon,"  Technical  Report  - 

ECOM  U.S.  Army  Electronics  Command,  Port  Monmouth,  NJ  (1978). 

60.  We  remark  that  these  N  vs  cp  curves  for  annealed  samples  sometimes 

ss  s 

contain  peaks,  e.g.,  see  Ref. 33.  We  treat  only  those  curves  without 
peaks  here  because  the  presence  of  the  peaks  causes  certain  difficulties 
with  the  expression  for  f®. 

61.  It  should  be  noted  that  a  term  A  (T)  could  appear  in  the  expression  for 
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f  in  (4.23).  However,  such  a  term  is  not  required  for  the  treatment  of 
any  static  or  dynamic  electrical  problems  involving  semiconductors. 

62.  The  particular  form  chosen  for  f  should  be  regarded  merely  as  a  good 
representation  of  some  ideal  function  truly  characteristic  of  surface 
response.  If  too  many  terms  in  the  power  series  in  cpg  had  been 
required,  different  (better)  functional  forms  with  fewer  defined 
material  surface  coefficients  would  have  been  considered. 


63.  In  obtaining  (4.25)  additional  assumptions,  i.e.,  Eqs. 0.24),  (4.24) 
and  the  particular  expression  used  for  (3.45)  in  '4.19),  were  employed. 
Therefore,  the  B^  coefficients  depend  on  parameters  independent  of  surface 

behavior  and  could  change  even  when  the  actual  surface  response  remains 
the  same.  In  contradistinction,  the  A^  are  fundamental  material  param¬ 
eters  descriptive  of  the  surface  response  alone. 

64.  b.E.  Deal,  M.  Sklar,  A.S.  Grove  and  E.H.  Snow,  "Characteristics  of  the 
Surface  State  Charge  (Q_  )  of  Thermally  Oxidized  Silicon,"  J.  Electrochem. 
Soc.,  114,  266  (1967). 
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65.  K.  Deuling,  E.  Klausman  and  A.  Goetzberger,  "Interface  States  in  Si-Si02 
Interfaces,"  Solid-State  Electron.,  15,  559  (1972),  Pig. 5. 

66.  In  calculating  Aq  from  the  measured  B  it  is  necessary  to  know  cpe 
and  cpox<  Values  for  these  quantities  may  be  found  in  Appendix  B. 

67.  And  hence  it  suffers  from  the  deficiency  noted  in  Ref. 63. 

68.  M.  Kuhn,  "A  Quasistatic  Technique  for  MOS-C-V  and  Surface  State 
Measurements,"  Solid  State  Electron.,  ^3,  873  (1970). 

69.  Electric  fields  of  this  magnitude  are  well  within  the  regime  of  field 
dependent  mobilities  in  silicon. 

70.  S.  Kar,  "Interface  Charge  Characteristics  of  MOS  Structures  with 

Different  Metals  on  Steam  Grown  Oxides,"  Solid  State  Electron.,  18, 

723  (1975),  Pigs. 5  and  Table  2.  This  is  a  case  where  data  is  given 

only  in  the  central  cp  region  and,  as  previously  mentioned,  the  field 
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term  is  well  represented  by  power  series  terms  in  cp  in  this  range. 
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As  a  result  the  coefficients  presented  here  are  composites  of  actual 

cp  terms  appearing  in  (4.23)  and  cp  terms  resulting  from  a  power  series 
s  s 

expansion  of  the  field  term. 

71.  The  present  macroscopic  theory  represents  a  thermodynamic  system  and 
hence  is  subject  to  the  laws  of  thermodynamics.  In  this  vein,  see  the 
introductory  sections  of  C.  Herring  and  M.H.  Nichols,  "Thermionic 
Emission,"  Rev.  Mod.  Rhys.,  _21,  185  (1949). 

72.  S.R.  DeGroot  and  P.  Mazur,  Non-Equilibrium  Thermodynamics,  (North 
Holland,  Amsterdam,  1962). 

73.  C.  Truesdell  and  R.A.  Toupin,  "The  Classical  Field  Theories,"  in 
Encyclopedia  of  Physics,  ed.  S.  Flligge  (Springer-Verlag,  Berlin,  I960), 
Vol.III/1,  Sec. 255. 

74.  The  question  arises  as  to  why  chemical  potentials  were  not  employed  from 
the  outset  here  since  much  of  the  existing  work  on  semiconductors  is 
grounded  either  explicitly  or  implicitly  in  the  chemical  potentials. 

The  reason  is  that  in  the  fully  macroscopic  balance  law  approach  used 
here  forces  and  pressures  are  fundamental.  Only  after  the  governing 
equations  are  obtained  are  the  pressures  abandoned  in  favor  of  more 
convenient  variables. 

75.  'rtiis  procedure  has  been  termed  the  principle  of  equipresence  by 
Truesdell.  See  Ref. 72,  Sec. 293. 
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76.  As  was  pointed  out  in  Ref. 55  detailed  specifications  of  cp  and  cp  are 
necessary  when  dynamic  problems  are  to  be  treated. 
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77.  The  expressions  in  (B.4)  can  be  found  in  Ref.l,  Ch.10.  It  is  noted 
that  if  we  were  to  choose  expressions  for  ee  to  match  those  in  the 
conventional  theory  for  degenerate  cases  then  cpe  would  be  taken  as 

-  (1/q)  [E  +  §(n/N  )]  (where  £  (n/N  )  is  defined  in  Appendix  II  of  Ref.l). 
c  c  c 

A  £ 

78.  Nhen  spatially  inhomogeneous  cases  are  considered  it  turns  out  that  a  cp 
different  from  cpe  can  be  defined  such  that  (l/pe)7pe  =  Vcpe.  This  cpe  in 
the  "ideal"  gas  case  is  related  to  cpe  by  ^  =  cpe  +  Ec/q.  In  the  event  that 
Ec  is  constant  everywhere  except  for  jumps  in  value  across  interfaces  - 
which  is  the  only  case  considered  here  -  both  cp®  and  cpe  satisfy  the  same 
differential  equation  (3.26)  on  either  side  of  apy  interface,  but  the 
boundary  condition  (3.40)  should  be  written  as  [cp®]  =  fe.  Since  in  this 
Appendix  we  are  interested  in  writing  this  boundary  condition  for  a  situ¬ 
ation  where  Ec  suffers  a  jump  in  value,  we  work  here  in  terms  of  q^.  Of 
course,  all  the  corresponding  expressions  for  holes  are  treated 
analogously. 
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79.  Ref. 5,  Ch.12.  Note  that  in  terms  of  cp  ,  (B.8)  may  be  written  cp  +  cp 

=  (-l/q)E^n\  i.e.,  (-l/q)E^n^  defined  here  is  equal  to  the  "electro¬ 
chemical"  potential  in  accordance  with  the  conventional  terminology. 

80.  Ref.l,  Ch. 10. 

81.  Ref. 7,  p. 468. 


Table  1  The  Material  Surface  Coefficients  Found  for  the 
Si-Si02  Interfaces  Studied  in  Ref. 70  where  Dif¬ 
ferent  Front  Contact  Metals  were  Used  on  Other¬ 
wise  Identical  Samples 


Contact  Metals 

Cr 

A1 

Mg 

(cm) 

-  4.  IX  10"5 

-  5.0 X  10"5 

-  7.0 X  10-5 

A2  (cm) 

1.4 X  10~4 

1.7 X  10-4 

2.4 X  10“4 

A3 

-  2.8 

-  2.5 

-  3.0 

A.  (volts  1) 

4 

-  2.1 

-  1.7 

-  2.8 

A_  (volts  2) 

D 

-  3.6 

-  2.9 

-  3.9 

FIGURE  CAPTIONS 


Figure  I  A  Schematic  Diagram  of  the  MOS  Structure 

Figure  2  A  Comparison  of  Experimental  Data  from  Ref. 64  with  the 

Theoretical  Curve  Obtained  Using  the  Expression  for  f® 

Given  in  Eq.  (4.23).  The  coefficients  are: 

B  =1.6X10  ®  coul/cm2  (assumed),  B  =+1.2X10  ^4  farad/cm, 

®  —8  2  ^  _7  2  2 
B_  =-4,3X10  farad/cm  ,  B_  =-1.4X10  coul/V  cm 

^  -7  3  2  J 

B^  =-5.2X10  coul/V  cm  .  The  A^  coefficients  are: 

A  =  -  7.6  X  10~5  cm,  A  =2.4X10~4  cm,  A,  --1.1X101, 

1  i  _i  1  2  -2  J 

A4  =-3.3X10  V  ,  A5  =-1.2X10  V  . 

Figure  3  A  Comparison  of  Two  Theore  tical  Curves  Obtained  Using 

Eq. (4.25)  with  Data  from  Ref. 68.  Curve  1  is  obtained  in¬ 
cluding  the  field  term  while  Curve  2  is  obtained  without 
the  field  term.  The  non-zero  B.  coefficients  are: 

—12  ^  —8  2 
Curve  1:  B.  =  1.9 X  10  farad/cm,  B  =  3.4 X  10  farad/cm  , 

-8  2  2  ^  -8  2 
B3  =  8.0X10  coul/V  cm  .  Curve  2:  B2=  -  6.4X10  farad/cm  , 

B3  =-4.4X10  7  coul/V2  cm2,  B4  =-1.8X10  ^  coul/V3  cm2, 

B&  =  -  3.1X  10"6  coul/V4  cm2 ,  Bg  =  -  2 . 1  X  10~^  coul/V5  cm2 . 

A  Comparison  of  Theoretical  Curves  with  Experimental  Data 
from  Ref. 70  Showing  the  Variation  of  Response  with  the  Use 
of  Different  Front  Contact  Metals  on  Otherwise  Identical 
Samples.  The  non-zero  material  surface  coefficients  for 
each  case  are  given  in  Table  1. 
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The  integral  form  of  the  conservation  of  energy  for  the  three -component 
combined  continuum,  consisting  of  the  two  conducting  fluids  plus  the  polarizable 
lattice,  states  that  the  time  rate  of  increase  of  total  stored  energy  in  the 
three  continua  in  any  fixed  volume  is  equal  to  the  rate  at  which  work  is  done  by 
the  fluid  pressure  forces  minus  the  rate  at  which  stored  energy  flows  out 
of  the  fixed  volume  plus  the  rate  at  which  energy  is  supplied  to  the  combined 
semiconducting  continuum  by  the  electric  field,  which  enables  us  to  write 


(e  +  pSee  +  p^e*1)  av 


I*  r  e  eh  h  eee 

[-p  n  •  v  -  p  n  •  v  -  n*v  p  e 

J  <V  <V  'V  A/  r 

s 

-  n  •  vhpheh]dS  +  j  EdV, 

V 


(A.  1 ) 


where  e  is  the  stored  internal  energy  per  unit  volume  in  the  lattice  con- 
.  e  h 

tmuum,  e  and  e  are  the  stored  internal  energies  per  unit  charge  of  the 

•  ♦  g  g  h  h 

conduction  electronic  and  hole  fluids,  respectively,  -p  n  •  v  and  -p  n  •  v 
are  the  rates  at  which  work  is  done  per  unit  area  by  the  pressures  acting 
on  the  two  charged  fluids,  n  *vepeee  and  n  •  vhpheh  are  the  rates  of  efflux 
of  stored  internal  energies  of  the  two  fluids  from  the  fixed  volume  and  E  is 
the  rate  of  supply  of  energy  per  unit  volume  to  the  combined  semiconducting 
continuum  from  the  electric  field.  The  quasistatic  Maxwell  electric  field 
acting  on  all  charged  elements  in  the  macroscopic  model  of  Section  2  supplies 
energy  according  to 

e  e  h  h 

E  =  p  E  •  v  +  p  E  •  v  +  E  •  —  ,  (A. 2) 

where  the  third  term  containing  the  electric  polarization  P  would  arise  from 

'V 

a  detailed  treatment  of  polarization  in  the  macroscopic  model^  but  is  here 
merely  taken  as  a  plausible  assumption. 


a 


A2. 


Taking  the  time  derivatives  in 

the  divergence  theorem,  (3.5),  (3.6), 

the  defined  material  derivatives 

_d_  _d_  d^_  e 

dt  =  at  »  dt  =  dt  ~  *  ~  * 


(A. 1),  substituting  from  (A. 2),  employing 
(3.17),  (3.18),  the  arbitrariness  of  V, 


£1 

dt 


_d 

dt 


=  rr  +  v  •  v 


(A.  3) 


and  in  accordance  with  the  model, making  the  plausible  assumptions 


e  e.  e  h  h  h 

e  =  e  (p  ,T)  ,  e  =  e  (p  ,T)  , 


(A.  4) 


we  obtain 


e  36®  L  h  deN  dT 
P  P  dT  /  dt 


■3 „  /  ~  _e  e  /  ,  ^  n  u*  n  * 

Ss.  +  (pe  _  p_\  d_2_  +  /  h  de_  _  £_\  a_g_  +  / 

dt  v  e  e/  dt  v  ^  h  h/  dt  v 
dp  p  dp  p 

_  ee  e  hh  h  e/pe  e\  h/p“  h\ 

’  5  at  =‘p  5,  -X  “P  £  -X  “Y  (“  +  e  huV  «  )  •  <A*5) 


BT 

h 


where  we  have  made  use  of  the  fact  that  VT  is  =  0  since  T  is  uniform  in  this 

•v 

work.  Equation  (A. 5)  is  commonly  called  the  first  law  of  thermodynamics  for 

the  thermodynamic  system  under  consideration. 

Since  the  first  law  of  thermodynamics  for  the  combined  continuum  is 

of  the  form  shown  in  (A. 5)  and  only  the  terms  on  the  right-hand  side  are 

dissipative,  the  expression  of  the  second  law  of  thermodynamics  may  be 
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written  in  the  form 


de 

dt 


+ 


3e 

dp 


e  de_ 
dT 


+  P 


h  deh\  dT 
dT  )  dt 


(p 


h  de_ 
dp* 
dP 
*  dt 


} 


(A.  6) 


where  T|  is  the  entropy  per  unit  volume.  Then, for  a  uniform  temperature 
state,  the  entropy  inequality  takes  the  form 


A3. 


^+ee)  +  Yh 

pe 

and  dT]/dt  is  the  rate  of  entropy  production. 

Equations  (A.6)  and  (A. 7) , which  arise  from  the  laws  of  thermodynamics 
as  applied  to  this  continuum,  are  the  basis  for  the  generation  of  two  types  of 
constitutive  relations s  those  from  the  former  are  termed  recoverable  and 
from  the  latter,  dissipative. 

In  obtaining  the  recoverable  constitutive  equations  it  turns  out  to 
be  convenient  to  define  a  thermodynamic  state  function  X  by, 


dTL 

dt 


-?[■ 


e_e 
p  E 


•  v®  +  phEh 

r  <v 


h  e/ 

X  +Y  { 


(4  ♦  <h)] 


2  o 


(A.  7) 


X*e  -  E  •  P  -  T|T , 


(A.  8) 


the  substitution  of  which  in  (A.6)  yields 

e  dee  =e  e  '  ^  '  h  -h'  h 


dE 


+  ( ne  2S-  -  d  p  +  (nh  _  Sl7\  d  p  +  P  . 

dt  +  Vp  7T  e)  at  +  Vp  7T  b;  at  + 1  at 


ap 


Sp 

e 


+  (-n  4.  ne  SSI  +  nh  ®  -  o 

+  \T1  +  P  ^F  +  P  Try  at  ~  0  • 


(A.  9) 


Since  (A. 9)  is  a  state  function  equation,  we  must  have 

X»X(E,T).  (A.  10) 


Taking  the  total  time  derivative  of  (A. 10)  and  substituting  it  into  (A. 9),  we 
obtain 


dE 


%jX  +  P> 


^  e  e.  ,e  e 

e  5e  _  p  \  d  p 

.  (  h  3eh 

hN 

-  2-) 

_h  h 
d  P 

.  e  ej  at 

ap  p 

+  \p  “T 

ap 

h) 

P 

dt 

,  -  ,  e  .  h 

♦n  +  p  -35- +  » 

deN  dT 

n 

Q> 

►3 

0.1 

rt-| 

1 

u  • 

(A. 11) 


Since  all  the  time  derivatives  appearing  in  (A. 11)  are  independent  and  (A. 11) 
holds  for  arbitrary  values  of  those  time  derivatives,  each  coefficient  must 


A4. 


vanish  separately,  which  yields 


5)1  e  5e6  h  3eh 
"  “  dT  “  p  3T  "  P  9T  ’ 


e  e,2  de  h 

P  =  (p  )  —r,  P  . 

ap 


,  h.  2  9eh 
(P  )  ~h 
3p 


(A.  12) 

(A.  13) 


Equations  (A. 12)  and  (A. 13)  are  the  recoverable  constitutive  relations  for 

6  1^1  € 

P,  Tj,  p  and  p  .  The  substitution  of  the  functional  dependences  for  x,  e  and 

e*1  shown  in. (A. 4)  and  (A. 10)  in  (A. 12) ^  and  (A. 13)  and  the  appropriate  use 

of  the  definition  of  D  give  rise  to  the  functional  forms  in  (3.19)  and  (3.21). 

<v 

As  usual,  chemical  potentials  are  defined  in  terns  of  e6  and  e*1  as^3 


e  _  3(oee6)  h  _  b(pheh) 

cp  =  -  *  V  ,  cp  , 

9p  3p 


(A.  14) 


0  II 

where  cp  and  9  are  the  chemical  potentials  of  the  conduction  electronic 

and  hole  fluids,  respectively.  These  are  clearly  related  to  the  fluid 
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pressures  through  (A. 13)  .  A  useful  relation  between  the  gradients  of  the 

chemical  potentials  and  of  the  pressures  in  the  respective  fluids  is  readily 
found  by  taking  the  gradients  of  (A. 13)  and  (A. 14),  with  the  results 


1  e  e  1  h  h 

—  ZP  -XfP  >  —  XP  -Sp  > 
p  p 


(A. 15) 


which  have  been  used  in  Section  3. 

Equations  (A. 13)  and  (A. 14)  enable  us  to  write  Eq. (A. 7)  in  the  form 

(A. 16) 


e  e  e  hh  h  ee  h  h 
p£’X+pS*X+y9+y<p  so 


Since  9®  is  a  function  only  of  pe  and  T,  and  cp*1  only  of  p*1  and  T  and  (A.  16) 
must  always  be  true  in  order  that  the  laws  of  thermodynamics  not  be  violated, 


we  must  have  the  dissipative  constitutive  relations 


A5. 


_e  _e.  e  h  e  h 

5  ml  (P  »P  >Z  >Z  >T)  » 

h  h ,  e  h  e  h 

E  -E  (p  ,p  ,v  v  ,T)  , 

e  e.  e  h  e  h 

-y  =  v  <p  ,p  ,  v  ,v  ,T)  , 

h  h .  e  h  e  h 

V  =  y  (p  ,P  ,  v  ,  v  ,T)  , 


(A. 17) 


along  with  certain  conditions  dictated  by  the  inequality  in  (A. 16)  which 

are  not  of  particular  interest  to  us  here. 

However,  since  the  recoverable  constitutive  equations  (A. 12)  and  (A. 13) 

depend  on  E,  there  is  no  logical  reason  to  exclude  E  from  the  dissipative 
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constitutive  equations  in  (A. 17)  .  Also,  since  the  fluids  have  been  chosen 

e  Ji 

to  interact  solely  through  the  generation/recombination  terms  y  and  y  ,  we 
exclude  the  dependence  of  ES  on  p*1  and  vh  and  of  Eh  on  p6  and  v6.  The  func- 
tional  dependence  of  the  dissipative  constitutive  equations  may  then  be 
written  in  the  form 

e  _e .  e  e  _  _h  _h ,  h  h  _ 

l  -£  (P  >Z  >  5  "S  (P  >£>T)  > 


e  e ,  e  h  e  h  _  h  h.  e  h  e  h 

y  =  Y  (P  ,p  ,v  ,v  ,E,T)  ,  Y  -y  (p  ,p  ,v  ,v  ,E,T)  , 


(A.  18) 


which  were  given  in  (3.20)  and  (3.22). 


A 


APPENDIX  B 


SOME  RELATIONS  BETWEEN  THE  MACROSCOPIC  DESCRIPTION 
AND  CONVENTIONAL  SEMICONDUCTOR  THEORY 

Ihe  macroscopic  theory  can  be  further  specified  (beyond  Section  4)  by  making 

0  J) 

particular  choices  for  the  energy  functions  e  and  e  which,  as  usual,  will 
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However,  since  ee  and  e*1  are  de- 


contain  macroscopic  material  parameters 
scriptive  of  the  conducting  fluids  in  the  interior  of  the  semiconductor  these 
parameters  may  be  related  to  and  calculated  with  reasonable  accuracy  from  the 
usual  quasi-microscopic  model  based  on  a  semi-classical  view  of  the  infinite 
periodic  lattice.  By  making  the  specific  selections 


(B.l) 

(B.2 


where  N  and  N  are  the  quasi-microscopically  defined  "effective  densities  of 
C  V 

states"  in  the  conduction  and  valence  bands,  respectively,  we  can  deduce  the 
Maxwell  gas  approximation  for  the  electrons  and  holes  in  the  semiconductor^. 
Equations  (B.l)  and  (B.2)  give  the  energies  per  unit  charge  for  this  particular 


model.  The  use  of  (B.l)  and  (B.2)  in  (A, 13)  and  (A. 14),  respectively,  allows  us 

h 


to  obtain77 


e  _  ^kTg 

q 

E 
_c 

q 


.  _  h  kTp  _ 
kTn  ,  p  =  -~r~  =  kTp  , 
'd 


(B,  3) 


**  =  -  log(_£L_)  =  _  _c  ±1  log(jL\ 

q  q  *\-qN  /  q  q  'Nc' 

.  s  *  -  i<4f) , 

q  q  \qN  /  q  q  \N  /  * 


(B.  4) 


where  n  and  p  are  the  electron  and  hole  number  densities,  respectively. 
Differentiating  (B.4)  with  respect  to  the  charge  densities  and  comparing 
the  result  with  (3.30),  we  can  obtain  the  Einstein  relations  for  the  Maxwell 


gas  model 


B2. 


When  no  electron  current  is  flowing  E6  =  0  and  hence  (3.26)  can  be 

78 

integrated,  with  the  aid  of  (B.4)^,  to  yield 

cp+®e-^  -  -p  lo^r)  >  <B*5> 

c 

where  n^  is  the  electron  number  density  in  the  bulk  of  the  semiconductor 
and  cp  has  been  assumed  zero  in  the  bulk.  In  the  conventional  quasi- 
microscopic  theory  for  nondegenerate  situations  at  thermal  equilibrium  an 
expression  for  n^  may  be  written  in  terms  of  the  Fermi  energy  E  ,  and  the 
conduction  band  energy  E^.  The  expression  is 

exp[(EF-Ec)/kT]  .  (B.6) 


Consequently,  when  no  electron  current  flows  we  may  combine  (B.5)  and  (B.6) 
to  obtain 

cp  +  cp®  *  |  (Jrc"V  ,  (B-7) 

where  it  is  to  be  noted  that  when  electron  current  flows  the  left-hand  side 
of  (B.7)  remains  well  defined  while  the  right-hand  side  does  not.  However, 
(B.7)  may  be  extended  to  include  current  flow  situations  by  employing  (B.7) 
to  define  a  quasi-Fermi  level79  for  electrons  E^11^ ,  with 


1 

q 


l 

q 


E 

c 


*e 

-  CP  -  cp  . 


(B.  8) 


The  relations  (B.7)  and  (B.8)  with  (B.4)^  substituted  are  in  agreement  with 

80 

expressions  in  the  literature  . 

Equation  (B.7)  [or  (B.8)]  may  be  used  to  exhibit  the  connection  between 
the  boundary  condition  (3.40)  and  the  conventional  semiconductor  nomenclature. 
For  the  semiconductor- vacuum  boundary  (an  oxide  interface  is  treated  ana¬ 
logously)  we  have  (for  the  no  current  situation), 


B3. 


rcpei 


«e  «e 

!  Cp  -cp 

^sera  Yvac 


-  (E  -  E) 
q  c  F 


■<P. 


(B.9) 


where  a  has  been  taken  to  be  identically  zero  and  cp  is  the  surface 
^vac  s 

potential.  We  may  rewrite  (B.9)  in  the  form, 


-  (E*  -  E  )  =  i  T-Xsr, 
q  c  P  q  ^  9 


(B.  10) 


where  E  is  the  "bent"  conduction  band  energy,  ¥  is  the  work  function  and 
c 

X  is  the  electron  affinity.  Equation  (B.10)  provides  a  ready  understanding 

e  . 

of  f  in  the  conventional  semiconductor  terminology.  When  an  n-type  semi- 

conductor  is  in  accumulation  T  is  small  and,  hence,  f  is  small.  This  reflects 

the  fact  that  in  accumulation  the  energy  of  the  near  surface  electrons  has 

been  lowered  and, hence, less  force  need  be  exerted  by  the  lattice  to  keep  those 

electrons  in  the  crystal.  In  inversion  the  reverse  is  true.  The  effects  of 

changes  in  temperature  may  be  argued  in  a  similar  manner.  It  is  important 

to  emphasize  that  while  (B.9)  or  (B. 10)  allow  the  interpretation  of  f 

within  the  framework  of  conventional  semiconduction  theory,  in  order  to  use 

(3.40)  in  the  solution  of  semiconductor  boundary  value  problems  it  is  neces- 

0 

sary  to  provide  an  expression  for  f  in  terms  of  the  dependent  variables  as 

is  done  in  Section  4.  For  holes  expressions  analogous  to  (B. 5)  -  (B.10) 

80 

may  easily  be  found  and  are  in  agreement  with  the  literature 

Wte  conclude  this  Appendix  by  utilizing  the  above  expressions  to  obtain 
relations  for  the  two  constants  cp®  and  $ox>  which  appear  in  (4.21)  and  are 
required  for  the  conversion  of  the  coefficients  in  (4.26)  to  the  funda¬ 
mental  coefficients  in  (4,23).  Since  no  current  flows  (B.7)  gives 
immediately 


cp 


J,  w 

q  (EC  -  V 


55  V  in  intrinsic  silicon  . 


(8.11) 


•■ii  a  *  «nfr'ii' 


B4. 


Concerning  the  oxide,  we  view  insulators  simply  as  semiconductors  with 
zero  mobilities.  As  a  result  of  this  assumption,  with  (3.23)^  and  (3.26) 

we  can  assume  uniformity  of  $e,  i.e. ,  cpe  =  =  const,  even  when  a  field 

,  ...  81 

is  applied.  Hence  , 

^  (E°X-E°X)  *  4.0  V  in  silicon  dioxide.  (B.  12) 
Expressions  for  and  nay  be  found  in  a  similar  manner. 


